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Galactic rotation curves and strange quark 
matter with observational constraints 

M. Kalam^,F.Rahaman^,Sk. M. Hossein^ and J.Naskar'^ 

Abstract 



We obtain the space time of the galactic core in the framework of general relativ- 
ity by taking the flat rotational curve as input and considering the matter content in 
the galactic core region as strange quark matter. We also obtain the energy density, 



radial and tangential pressures of the core strange quark matter. Significantly we 
have shown that Bag Constant takes an important role to stabilize the circular orbit 
of the test particles. We also give a limit of the Bag Constant for the existence of 
• quark matter in the galactic halo region. 
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Presence of dark matter was suspected first by Oort^ [21 E]- Later Zwicky [U E] con- 
cluded that a large amount of non-luminous matter remains hidden in the galactic haloes. 
Observation of flatness of galactic rotation curves [6l [71 [H [9l [TOl [11] was consequitively 
confirmed the Zwicky's suggestion regarding the presence of dark matter. Afterwards, 
gravitational lensing of objects like bullet clusters and the temperature distribution of 
hot gas in galaxies and galactic clusters have further confirmed the existence of dark 
matter. The flatness of the galactic rotation curves indicates the presence of much more 
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matter within the galaxy than the visible material. Estimation supports a spherical dis- 
tribution of matter (galaxy) is surrounded by a dark matter halo. Gravitational effect of 
dark matter is more manifest at larger radius[T2]. Now a days, the amount of dark matter 
present in the Universe has become known in more precise manner. Cosmic Microwave 
Background Radiation data indicates that nearly 85% of total matter in the galaxy is 
dark in nature. Big Bang analysis and cosmological observational data indicates that the 
bulk of dark matter may be cold or warm,stable or long-lived and non-interacting with 
the visible matter. However, after a long and through investigations, a little is yet known 
about the nature of dark matter. 

Gell-Mann jl3 j and Zweig [14] had independently proposed that hadrons are consists of 
even more fundamental particles called quarks, a proposition which got experimental sup- 
port later on. As quarks are not free particles, the quark confinement mechanism have 
been dealt with great details in Quantum Chromo Dynamics (Q CD). In the MIT bag 
model[T5],it was suggested that the quark confinement is due to a universal pressure Bg, 
which is called the bag constant. Farhi and Jafee |T6] and Alcock et al [IT] had shown 
that the value of the bag constant, Bg should lie between 60-100 MeV / frn? , for a strange 
quark matter. 

There are a number of proposed candidates for dark matter. One of them is standard 
cold dark matter(SCDM)p^ 119]. In this hterature, we suggested that quark matter may 
be such candidate which is also proposed by several researchers [201 EH 122] • Also, it is 
well accepted that quark matter is exist at the centre of neutron stars, strange stars [201 
[23| [211 [251 [23 [23 [21] • immediately after the Big Bang the Universe underwent a Quark 
-Qluon-Plasma(QGP) phase. In the famous Large Hadronic Collider (LHC) experiment. 
Scientists recreate the conditions similar to those encountered before and in the early 
hadronization period[2S].As the expansion takes place, the Universe cools down and the 
hot Quark-Gluon-Plasma(QGP) freezes slowly to produce individual hadrons [50] . Most 
acceptable theory of strong interaction. Quantum Chromo Dynamics (QCD), suggests that 
under extreme condition a hadronic system can undergoes a phase transition from confined 
hadronic matter to the QGP phase. It is suggested that the core of neutron stars may 
consists of cold QGP[31J. 

The purpose of the present work is to show that, with the input of flat rotation curve and 
assuming the dark matter contents as strange quark matter, the gravity in the halo region 
is attractive under certain condition where Bag Constant takes an important role. We also 
want to investigate the presence of quark matter in the galactic halo region. This is the 
main motto of our work. 



2 



2 Field equations and general results 



The general static spherically symmetric spacetime is represented as 

ds^ = -e'^'Ut^ + e^^'Ur^ + r\de^ + sin^ed(P^) (1) 

where iy{r) and A(r) are the metric potentials and are function of the space coordinate r 
only. We further assume that the energy-momentum tensor for the strange quark matter 
filling the interior of the galaxy may be expressed in the standard form as 

Tij = diag{p, -pr, -pt, -Pt) 

where p ,pr and pt correspond to the energy density, radial pressure and transverse pressure 
of the baryonic matter, respectively. 

The Einstein's field equations for the metric (1) in presence of strange quark matter are 
then obtained as (with G = c = 1 under geometrized relativistic units) 

e-" (--\)+\ = 8np, (2) 

(7 + ^) - = ^^Pr. (3) 

-e-(." + -- — + ^j=8.,. (4) 

Following the MIT bag model, we express the strange quark matter equation of state(EOS) 
in the form 

Pr = \{p-^B,), (5) 

where Bg is the Bag constant (in MeV/ fm? units). 



3 Solutions for the strange quark matter model 

Assuming the known flat rotation curve condition t><^ = [-^^^Y^"^ = constant tangential 
velocity |33l [32] obtained a solution of it as [31] 

= Bot\ (6) 

where / = 2f^ and Bq is an integration constant. 

It is to be noted that the observed rotational velocity v"^ becomes more or less a constant 
with v'l' ~ 10~^(300 km/s) for a typical galaxy [521 EH]- 
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Now, from equations (2)- (5) and then using the simple expression in the equation(6) we 
get the simphfied form as 



(e~^A') + (4 + 3/)^ = ^ - 327iBgr 



(7) 



Substituting e = x, we get 



, ax 4 
X H = — h cr 



(8) 



where a = 4 + 3/,c = —327iBg. Solving the above equation we get the metric potential as 



4a + car^ + 8 D 
a{a + 2) ^ 



where D is integration constant. 



(9) 



Therefore, with the fiat rotational curve condition, the metric (1) becomes 



4a + car^ + 8 D 
a(a + 2) r'* 



H -1 



(10) 



Using equation (6) and (9) in (2)-(5), we get 



1 (D{a - 1) 3c ^ a - 4 



57r 



a + 2 



Pr 



4a + 8 + car^ D 



a(a + 2) 



Pt = ^ 



2c 



aD 



4a + 8 + car2 D\ (1 1 

a(a + 2) ^ ^ V 4 ^ 2/ Va + 2 



(11) 
(12) 

(13) 



The pressure anisotropy is a good feature from the point of view of exterior matching. 
The solution can be matched to the Schwarzschild exterior metric at the boundary of the 
halo[38]- As the strange matters are not exotic in nature, they must satisfies the Null 
Energy Condition(NEC) [39]. We find that the NEC is satisfies at everywhere i.e. p > 
and p+Pr > 0. However, we note that here the Strong Energy Condition(SEC) is violated 
( see figure 1). Therefore, we comment the quark matter in the galactic halo region is 
marginally real like. 
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Figure 1: {Left) Variation of p with distance r(Mpc) for the specific values of the 
parameters, D = 10^'' , / = 10^^ , Bg = .00013. (Mz(i(i/e)Variation of p+Pr with distance 
r(Mpc) for the same specific values of the parameters. Variation of p + Pr + '2^Pt 

with distance r(Mpc) for the same specific values of the parameters. 

4 Few Features 



4.1 Stability 

Let, the four velocity, = ^ for a test particle moving solely in the galactic halo region 
(restricting ourself to 6 = 7i/2), the equation g^^U^U" = —m^ can be cast in a Newtonian 
form as 1371 



dr 
dr 



+ V{r) 



(14) 



which gives 
V{r) = 



-I \ ia+S+car^ , D l 
' ^ a(a+2) 'i 



Bo 



+ 



4a + 8 + car^ D 



a{a + 2) r" 



+ - 1 + ^ 



mo mo 



(15) 



(16) 



where the constants E and L, respectively, represents the conserved relativistic energy 
and angular momentum per unit mass of the test particle. Circular orbits of the test 



particles are defined by r = i? =constant, so that 
two conditions we get the conserved parameters as: 

L = ±^[^R 



m 

dr 



and f\r=R 



2-1 



0. From these 



17) 
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and using this in V{R) — —E'^, one can easily get, 



2-1 



(18) 



For stable orbit 



for L and E in 



'V\ 



—B. < and for unstable orbit 
we get the final results as 



d^V\ 



=R > 0. Putting the expressions 



Thus 



(py 



=R 



2/_ 



4a + 8 + cai?^ 
a(a + 2) 



+ 



+ 



2(2/ -I) [ aD 



i?V i?(2-/) V^"+^ « + 2 



2ci? 



(19) 



=R < so that circular orbits are stable under consideration when D > and 



rj . (a+2)(2-0 
^ 647r(2-2/-i2) 



_8Z_ I 



2D 



a(l-20 
2-i 



The first condition is always satisfied since, 

from equation(lO), we see that D has the same dimension as radius, R. Again, as the 
Bag Const., Bg and radius R are always positive , so D is obviously positive. As this 
results is in agreement with the observations, therefore we can conclude that the value of 



SI 



I 2D 



l- 



a(l-2l) 



2-1 



, if indeed. 



the Bag constant, should be less than 64^12^-?/-?^) 
there exists any quark matter present in the galactic haloes. Therefore from this model, 
we can give the limit and a rough estimation of the Bag constant within the galactic halo 
region. Thus we note that Bag Const., plays an important role to stable the circular 
orbit. 



4.2 Attraction in Strange quark matter 



Observational indication is that the gravity on the galactic scale is attractive. Again, the 
stable circular orbit indicates that the particles are being accelerated towards the Centre 
of the Galaxy. We can sec it by studying the geodesic for a test particle that has been 
placed at a circular path of radius r. 
Now, 

+ — ^ = 0. (20) 



dT dr 



This equation implies that 



dr2 



4a+8+car^ , _D 
o(a+2) 



Da 



2car 
a{a+2) 



Aa+S+car^ 
a(a+2) 



(21) 



Obviously, the quantities in the square brackets are positive. Thus particles are attracted 
towards the centre with some boundaries. 
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4.3 Gravitational Energy 



We can easily determine the total gravitational energy Eg between two fixed radii, say, ri 
and r2 



Eg = M - Em = 4:n I [1 - V^]pr^dr 
= Air 



r2 



4a+car2+8 , D 



a{a+2) 



+ — 



1 fD{a-l) 3c a- A 



where 



T2 



M = Att pr dr 



r'^dr. (22) 



(23) 



is the Newtonian mass given by 



M = 47r 



>ri 



T2 

pr^dr = - 
^ 2 



Da 



D 



(2 - a)r'^-2 (1 - a)r'^-i 



a - 4 
_l f 



cr 



r2 



(24) 



Here, Em is the sum of the other forms of energy like the rest energy, kinetic energy, and 
internal energy. 

For the specific values of the parameters ri = 1 , r2 = 3 , D = 10"'' , I = 10"^ , Bg = 
.00013, we calculate the numerical value of the integrand (22) describing the total gravi- 
tational energy as Eg = —0.00168098, which indicates that Eg < 0, in other words, there 
is an attractive effect in the halo region. This result is very much expected for the matter 
source in the galactic halo region that produces stable the circular orbit. 



4.4 Observational Constraints 



We can rewrite the metric(l) in the form 



ds^ = -e^^^^^dt^ + 



2m{r) 



dr' + r\de' + sin'ed<i)') 



(25) 



where, the metric functions are given by 



m{r) 



$(r) = - (InEo + /Inr) 

Aa + car"^ + 8 D 



1 - 



a(a + 2) 



(26) 
(27) 
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These two functions are not necessarily the same as the potential and mass functions are 
found from the observations. The first post-Newtonian approximation |10], the gravita- 
tional potential $(r), in general relativity, is given by 

V'* ^ Rtt ^ 47r(p + pr + 2pt). (28) 

which reduces to the equation 

V^$iv = 47rp (29) 

provided the contributions of the pressures are negligible in comparison to energy den- 
sity.Here, $7v indicates the Newtonian potential. This general relativistic $ also appears 
in the wavelength shifts z±. of an emission line of a massive particle for an edge-on galaxy 
is [m, 112] is given by 



1 



^1 -r<l''( 



± \ b \ ^/r^'(r 



(30) 



where ± signs indicates the approaching and receding particles and b is the impact pa- 
rameter. The above equation(30) can approximately written as |43] Zj^ ~ r$'(r). 
Therefore, the usual methods for obtaining potential for rotation curve measurements 
yields the pseudo-potential 

$RC = $ 7^ ^N. (31) 

From equation(28) one can say that pseudo-mass, 

f Ir 
rriRc = r^$'(r) ^ An / (p + Pr + 2pt)r^dr = — (32) 



The pseudo-mass muc becomes Newtonian mass M(r) when pressure contributions are 
negligible in comparison to energy density. 



We can also treat Photons as probe particles as they interacts with the gravitational field 
in the core of the galaxy during their travel to the observer. The gravitational effect on 
photon motion can be measured in terms of a refractive index /i(r)^44j.The geodesic of 
both massive and photonlike particles could be expressed exactly in terms of a single 

2 

generalized refractive index N = ^ where v is the three-velocity of the particle. For 
light motion, v = -. 

For unspecified metric functions $(r) and m(r),[l3] argue that /i(r) = 1 — 2$;ens + [$fens] 
where they defined the lensing pseudo-potential as 

1 , , If m(r) , , , 

^lens = ^Hr) + -J -j^dr. (33) 

Another pseudo-mass obtained from lensing measurement has been written as |13] 

miens = ^r2$'(r) + ^m(r). (34) 
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The first order approximation of Einstein's equation yield 

Pir) ^ ■^^[2ml^,{r) - m'jic{r)] (35) 
AnPipr + 2pt) ^ 2[m'j,c{r) - m;,„,(r)] (36) 

wliere tlie riglit-liand sides represents, respectively, pseudo-density and pseudo-pressures. 
Furthermore, [43j defined a dimensionless quantity 

where the pseudo-quantities on the right hand side of the equations (35-37) determine, 
respectively, the observed density, pressure and equation of state. 



It is to be mentioned that the observable quantities are the pseudo-quantities. The general 
procedure to determine the metric are as follows: Once anyone is able to observationally 
determine the profiles of pseudo-quantities, one can work backwards to find the corre- 
sponding metric functions. This kind of reverse technique observational astrophysicists 
use. If the observed pseudo-profiles fit (up to experimental error) with the analytic profiles 
of a priori given metric functions, one can say that the solution is physically sustainable. 
Otherwise, it has to be ruled out as non- viable. In this sense, the observed pseudo-profiles 
play the role of constrains on the possible metric solutions. 

For the presence situations, we obtain the following constraint equations on $(r),m(r), 
the pressure profile and the equation of state: 

$RC = ^(lnEo + /lnr) 



Ir 
~2 



lens 



[l + l 



4 , In r 1 „ 

-) h -InEn H 

4 4 ° 4a 



D 



ac 



2(a + 2)_ 



2(m', 



RC 



m 



lens) 



co{r) 



2im\ 



RC 



m 



mie 
I 

2 ~ 

lens/ 



+ 



r 



1 - 



4a 8 -f car^ D 



4 V' a(a + 2) 
Aa + d, + car'^ D 
a(a + 2) r'^ 



Aa+8+car^ 
a{a+2) 



D_ 



2 car 
a{a + 2) 

2car 
' a{a+2) 



Da 



+ 



Da 



a+1 



3 (2m; 



lens 



RC) 



4a+8+car^ 
a{a+2) 



D 



+ 



2car 



a(a+2) 



+ 



Da 



(38) 
(39) 
(40) 

(41) 

(42) 

3/ 
2" 

(43) 



It is to be noted that the pseudoquantities given in Eqs. (38)-(43) are actual observables 
from the combined measurements of rotation curves and gravitational lensing. Equation 
(43) provides a convenient parameter that gives a 'measure' of the equation of state of the 
halo field, calculated from a combination of rotation curves and lensing measurements. 
Figure 2 shows that the value of u is negative indicating repulsion. 
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Figure 2: The variation uj{r) with respect to r. 

5 Conclusions 

We obtain the space time of the galactic core in the framework of general relativity by 
assuming the fiat rotational curve condition and considering the matter content in the 
galactic core region as strange quark matter. The spacetime metric is neither asymtotically 
fiat nor a spacetime due to a centrally symmetric black hole. We also obtain the energy 
density, radial and tangential pressures of the core strange quark matter. We see that 
the NEC is satisfies. However, we note that the SEC is violated. Therefore, it is obvious 
that quark matter is not purely real like and it behaves like dark matter. Significantly 
we shown that Bag Constant takes an important role to stabilize the circular orbit of the 
test particles. We also give a limit of the Bag Constant for the existence of quark matter 
in the galactic halo region. 



Acknowledgments 



We are thankful to Dr. G.C.Shit for helpful discussion. MK and PR gratefully acknowl- 
edge support from the Inter-University Centre for Astronomy and Astrophysics(IUCAA), 
Pune, India under Visiting Associateship programme under which a part of this work was 
carried out.SMH is also thankful to lUCAA for giving him an opportunity to visit lUCAA 
where a part of this work was carried out. PR is personally thankful to PURSE, DST and 
UGC for providing financial support. 



10 



References 

[I] J. Oort, Bull. Astron. Ins. Nether V, 189 (1930). 
[2] J. Oort, Bull. Astron. Ins. Nether. V, 192 (1930). 
[3] J. Oort, Bull. Astron. Ins. Nether. V, 239 (1930). 
[4] F. Zwicky F., Helvet. Phys. Acta 6, 110 (1933). 
[5] F. Zwicky, ApJ, 86, 217 (1937). 

[6] K. C. Freeman, ApJ 160, 881 (1970). 

[7] M. S. Roberts and A. H. Rots A. H., A&A 26, 483 (1973). 

[8] P. Ostriker, P. J. E. Peebles, and A. Yahill, ApJ 193, LI (1974). 

[9] J. Einasto, A. Kaasik, and E. Saar, 1974, Nat 250, 309 (1974). 

[10] V. C. Rubin, N. Thonnard, and W. K. Ford Jr., ApJ 225, L107 (1978). 

[II] Y. Sofue and V. Rubin, ARA&A 39, 137 (2001). 
[12] K.G. Begeman, A&A 223, 47 (1989). 

[13] M.GellMann, Phys. Lett. 8, 214 (1964). 

[14] G. Zweig, Ed. Developments in the quark theory of hadrons. Lichtenberg 1, 22 (1964). 

[15] A. Chodos, R.L. Jaffe,K.Johnson,C.B.Thorn and V.F.Weisskopf, Phys. Rev. D9, 3471 
(1974). 

[16] E. Farhi and R. L. Jaffe, Phys. Rev. D30, 2379 (1984). 

[17] C. Alcock, E. Farhi and A. Olinto, Astrophys. J.310, 261 (1986). 

[18] G. Efstathiou, W. Sutherland, and S. J. Madox, Nat 348, 705 (1990). 
[19] A. C. Pope, et al., ApJ 607, 655 (2004). 

[20] M.A. Perez-Garcia,J.Silk,J.R.Stone,Phys. Rev. Lett. 105, 141101 (2010). 
[21] M. M. Forbes, K.Lawson, A.R.Zhitnitsky, Phys. Rev. D 82, 083510 (2010). 
[22] S. Weissenborn et al, Astrophys.J. 740, 114 (2011). 
[23] J.J.Drake, Astrophys.J. 572, 996 (2002). 

[24] P. Jaikumar, S. Reddy, A.W. Steiner, Phys. Rev.Lett. 96, 041101 (2006). 

[25] H. Rodrigues, S.Barbosa Duarte, J.C.T. de Oliveira, Astrophys. J. 730, 31 (2011). 

11 



[26 
[27; 
[28 
[29 
[30 
[31 
[32 

[33 
[34; 

[35 
[36 

[37; 

[38 
[39 
[40 

[41 
[42 
[43 
[44; 
[45 



G.H. Bordbar, A.R.Peivand, arXiv: 1103.1250. 
J. Madsen, Lect. Notes Phys. 516, 162 (1999). 
F. Weber, Prog. Part. Nucl. Phys. 54, 193 (2005). 

EMUOl Collaboration, M.I.Adamovich et al, Phys. Lett. B 263, 539 (1991). 
M. Fromerth, J. Rafelski, in: CINPP Conference, Kolkata, February 2005, p. 22. 
T. Walker et al, Astrophys. J. 376, 51 (1991). 

F. Rahaman, R. Mondal, M. Kalam, and B. Raychaudhuri, Mod. Phys. Lett. A 22, 
971 (2007). 

S. Chandrasekhar, Mathematical Theory of Black Holes (Oxford University Press). 

F. Rahaman, M. Kalam, A. DeBenedictis,A. A. Usmani, and S. Ray, Mon. Not. Roy. 
Astron. Soc. 389, 27 (2008); F Rahaman et al, Phys.Lett.B 694, 10 (2010). 

T. Matos, F. S. Guzman, and D. Nunez, Phys. Rev. D 62, 061301 (2000). 

P. Salucci and M. Persic, arXiv: astro-ph/9703027] 



K. K. Nandi, A. I. Filippov, F. Rahaman, S. Ray, A. A. Usmani, M. Kalam, and A. 
DeBenedictis, Mon. Not. Roy. Astron. Soc. 399, 2079 (2009). 

S. Bharadwaj and S. Kar, Phys. Rev. D 68, 023516 (2003). 

D. Hochberg and M. Visser, Phys. Rev. Lett. 81, 746 (1998). 

C. W. Misner, K.S. Thorne, J. A. Wheeler, Gravitation. Freeman, SanFrancisco,p. 601 
(1973). 

U. Nucamendi, M. Salgado, D. Sudarsky, Phys. Rev. D 63, 125016 (2001). 
K. Lake, Phys. Rev. Lett. 92, 051101 (2004). 

T. Faber and M. Visser, Mon. Not. Roy. Astron. Soc. 372, 136 (2006). 
P. Boonserm et al Class. Quant. Grav. 22, 1905 (2005). 

D. Lynden-Bell, J. Katz and J. Bicak, Phys. Rev. D 75, 024040 (2007). 



12 



